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Abstract 

Let X„ = (xi, . . . , x n ) and / £ R[X„, k]. The problem of finding all k such that f(X„, k) > for 
all X n £ R n is considered in this paper, which obviously takes as a special case the problem of 
computing the global infimum or proving the semi-definiteness of a polynomial. For solving the 
problems we propose a simplified Brown's CAD projection operator, called Nproj, of which the 
projection scale is always no larger than that of Brown's. For many problems, especially when 
n > 3, the scale of Nproj is much smaller than that of Brown's. As a result, the lifting phase 
is also simplified. Some new algorithms based on Nproj for solving those problems are designed 
and proved to be correct. Comparison to some existing tools on some examples is reported to 
illustrate the effectiveness of our new algorithms. 
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Introduction 



Let's consider the following three well-known problems. 

Problem 1. For X n G R" and / G R[X n ], prove or disprove VA„ G R n (/(X„) > 0). 
Problem 2. For X n G R" and / G R[X n ], find the global infimum inf /(R n ). 
Problem 3. For X n G R™ and / G R[X ni k], find all k G R such that VA„ G 
l n (f(X n ,k)>0). 

A lot of work has been done for Problem 1 since Hilbert (Hilbert 18881). For related 



classical results, see for example, Bernstein (1915); Artin (1927); PolyaJ (|L928j); Hardy 
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Problem 2 can be regarded as a generalization of Problem 1. Various methods based on 
different principles have been proposed for solving Problem 2, including methods based 



on Grobner base (Hagglof et al. 1995 



based methods ( Parrilo 2000 Lasserre 



Hanzon 2003), semi-definite programming, SOS 



2001 Jibetean et al. 2005 Nie et al. 2006 Safey 



2008| |Ha et al.[ |2009| |Guo et al.[ [2010J), and methods based on Wu's method ( |Xiao et 
al. 2011). Semi-definite programming is based on symbolic- numerical computation and 
it returns numerical solution. Furthermore, in some cases it may return value larger than 
the supremum. Some methods need additional assumptions, for example: assuming that 



the polynomial can attain the infimum ( Nie et al. 2006 ) or assuming that the zero of 



1995). 



the set of the first partial derivatives is zero-dimensional ( Hagglof et al. 

Problem 3 is more general. It is a typical problem of quantifier elimination on real 
closed field. Therefore, theoretically it is feasible to apply the algorithm of Cylindrical 



Algebraic Decomposition^ AD) (Collins 1975 1998 Collins et al. 1991). The original 



CAD algorithm is of high complexity since the algorithm process of CAD projection phase 
involves a large amount of resultant calculation and the lifting phase needs to choose a 
sample point in every cell. Hence a lot of work tries to improve the CAD projection. A well 
known improvement is Hong's projection operator which is applicable in all cases (iHong 



1990 ). For many problems, a smaller projection operator given by McCallum ( 1988 1998 1, 
with an improvement by Brown (2001 1, is more efficient. Strzebonski (2000) proposed the 
algorithm of Generic Cylindrical Algebraic Decomposition(GCAD) for solving systems 
of strict polynomial inequalities, which made use of the so-called generic projection, the 



same projection operator as Brown (2001) proposed later. Based on Wu's principle of 
finite kernel (Wu 1998 2003), Yang proposed without proof the successive resultant 
method (Yang 2001 Yang et al. 20081 to solve the global optimization problem, in 



which Brown's projection is used in the projection phase and only sample points from 



the highest dimensional cells need to be chosen in the lifting phase. McCallum (1993) 
once pointed out that in order to prove a polynomial inequality, we need to choose sample 
points only from the highest dimensional cells, but he did not use the Brown projection. 



Xiao (20091 proved that, in terms of the Brown projection, at least one sample point can 
be taken from every highest dimensional cell via the Open CAD lifting. The algorithms 
proposed in this paper will also use this conclusion. 

In this paper, we simplify the Brown projection in order to solve Problems 1, 2 and 
3 more efficiently. In fact, we propose a simplified Brown's CAD projection operator, 
called Nproj , of which the projection scale is always no larger than that of Brown's. For 
many problems, especially when n > 3, the scale of Nproj is much smaller than that of 
Brown's. As a result, the lifting phase is also simplified. Some new algorithms based on 
Nproj are designed and proved to be correct. Comparison to some existing tools on some 
examples is reported to illustrate the effectiveness of our new algorithms. In order to 
prove the correctness of our algorithms, we find an interesting property: for h, g G M[A ?l ] 
with gcd(h(X n ), g(X n )) = 1 and a connected open set U in K™, C/\Zero(/i, (?) is also a 
connected open set, where Zero(h,g) means the set of common real roots of h and g. 

The structure of this paper is as follows. Section 2 introduces basic definitions, lemmas 
and concepts of CAD and Brown's projection. Section 3 proves the correctness of the 



successive resultant method proposed in Yang (2001 ). In Section 4 and Section 5, our new 



2 



projection operator Nproj is introduced and a new complete algorithm based on Nproj 
is proposed for solving the above three problems. The correctness of our algorithm is 
proved. The last section includes several examples which demonstrate the process and 
effectiveness of our algorithm. 

2. Preliminaries 



In this paper, if not specified, for a positive integer n, X n , Y n , A n , B n and 0„ denote the 
point (x 1 ,x 2 ,...,x n ) £ R n , (yi,y 2 ,...,y n ) & K™, (<n, a 2 , • • • , a n ) e M", (61, b 2 , . . . , b n ) e 
R n , and (0, 0, . . . , 0) G R n , respectively. 

Definition 1. For X ni Y n £ W 1 , the Euclidean distance of X n and Y n is defined by 



Definition 2. For X n e R n , let Bx n (r) be the open ball which centered in X n with 
radius r, that is 



Definition 3. Let / € R[X„], the set of real zeros of / is denoted by Zero(/). Let L be 
a subset of K[X„]. Define 



The elements of Zero(L) are the common real zeros of L. 

Definition 4. The level for / <G M^n] is the largest j such that deg(/, Xj) > where 
deg(/, Xj) is the degree of / with respect to Xj. The level for / e M.[X n , k] is the largest 
j such that deg(/, Xj) > and the level is zero if / is a polynomial in k. 

Definition 5. For polynomial set L C M[X„] or L C R[X n , k], L % is the set of polynomials 
in L with level i. 



Definition 6. For X n ,Y n e W 1 , we denote by X n Y n the segment X n -4- Y n . For 
X%,X*,...,X%, we denote hyX° n -> X \ ^ the broken line through X°,X^,...,X, 



The following two lemmas are well-known results. 

Lemma 7. Let f,g € K[^n]; f/iene exist A, B e M[X n ] smc/i that Af + Bg = Res(f,g,x n ) 
with deg(A, x n ) < deg(g 7 x n ) and deg(B,x n ) < deg(f,x n ), where Res(f,g,x n ) is the 
resultant of f and g with respect to x n . 

Lemma 8. Let f(X n ) e K[X„] and r be a real positive number. If f = for all X n e 
B 0n (r), then f = 0. 

For convenience, we denote by (f,g) the greatest common divisor of / and g in R[X„]. 




Bx n (r) ■= {Y n G K" I p(X„,F„) < r}. 



Zero(L) = {X n e M n |V/ G L, /(X n ) = 0}. 



in turn. 
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Lemma 9. For f,g € WL[X n ], if f and g are coprime, then after any linear invertible 
transform, f and g are still coprime, namely for A e GL n (R), X* T = AX% + B T , then 
(f(X*),g(X*)) = l. 

Proof. If (f(X*),g(X*)) = h(X n ) and h is not a constant. Then h{A~ l {X n - B%)) is 
a non-trivial common divisor of f and g in M[X n ], which is a contradiction. □ 



Lemma 10. Suppose f,g£ R[X n ] and (/, g) = 1. For any X n _i e 



id r > 0, the 



exists X' n _ 1 e 



such that p(X n _ 1 , X^^) < r and for all x' n e R, (X^_ 1; a;^) ^ Zero 



(/,<?)• 



Proof. Otherwise, there exist X®_ 1 = (x^x®, ■ ■ ■ >#n_i) € R™ -1 and ro > 0, such that 
for any X n _\ = (x\, x<i, ■ ■ ■ , x n -i) satisfying p(X^_ 1 , X n _i) < ro, there exists an x n € K 
such that f(X n ) = g(X n ) = 0. Thus Res(f,g,x n ) = holds on B ro (X^_ 1 ). From Lemma 
8, we get that Res(f,g, x n ) — holds for all X„_i e M™ -1 meaning (f,g) is non-trivial, 
which is impossible. □ 



Definition 11. For / <G R[x], say 

f{x) = aix 1 H h a ,ai ^ 0. 

The discriminant of f{x) is 

discrim(/, x) = af~ 2 (zi - zjf, 

i<j 

where Zi (i — 1,2, ... ,1) are the complex roots of the equation f(x) = 0. 

The following well-known equation shows the relationship between discrim(/, x) and 
R*®(f,fL,x), 

a;discrim(/, x) = (— 1) 2 Kes(f,f x ,x). 
For real parametric polynomial f(X n ,x), if the leading coefficient lc(f,x) ^ 0, the dis- 
criminant of f(X n ,x) can be written as 



i(i-i) 

discrim(/, x) := (— 1) 



/l 


ai-i 




ai-2 
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If a; = a;_i = at point X„, from the above expression, discrim(/, x) equals at this 
point. 
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Lemma 12. (Weiss, 1963) Let fix) E R[x] be a monic squarefree polynomial of degree 
I, the sign of its discriminant is (— 1)t~, where r is the number of its real roots. 

It is clear that the conclusion of the above lemma still holds when lc(/, x) is positive. 

Lemma 13. Given a polynomial f(X n ,x n +i) E R[X n , x n+ i], say 

i 

f(X n ,x n+1 ) = ^2a. t xl l+1 ,ai ^ 0, 



where a^ is a polynomial in X n (i = 0,1, ... ,1). Let U be an open set in WL n . If f(X n , x n+ i) > 
on U xR, then there exists a positive integer m such that I — 2m and 

(— l) m discrim(/, x„+i) > and lc(f,x n+ i) > for all X n E U. 



Proof. Obviously, I = deg(/, x n+ i) is even and lc(/, x n +i) is positive semi-definite on 
[/xl. Let I = 2m for some positive integer m. For any given X n E U, if f(X n ,x n+ i) has 
square factors, then discrim(/, x n +\) — 0. Since / is positive semi-definite, if a; = at 
X n , then a;_! = at X n . Thus discrim(/, x n+1 ) = at X n . Now, suppose / is squarefree 
and ai ^ at X n . Because / is positive semi-definite, / has no real roots. By Lemma 
12, the sign of discrim(/, x n+ i) is (— l) m at X n . That completes the proof. □ 



Definition 14. Suppose hE M.[X n ] can be factorized in R[A n ] as: 

h = hi ll h 2 r ' 1 . . . h m tm , 
where hi(i = 1, . . . , m) are pairwise different irreducible polynomials in 

m 

sqrfree(ft.) = hi. 



n ]. Define 



Lemma 15. (Yang et all 2008) Given a real polynomial f with real parameters, say 



/(a, a;) 



-\x 



where el — [a m , ...,oq) is real parameter. Let i?(a) = sqrfree(Res(/, /', x)). If c\ and 
ci are two points in the same connected component of parameter space R(a) =/= 0, then 
/(ci,x) and f(c%,x) have the same number of real roots yi(ci) < 1/2(01) < . . . < yd{ci) 
and yi (C2) < 2/2(02) < • • ■ < Dd{ c 'i)- Moreover, j/j(a)(i = 1, 2, . . . , d) is continuous in the 
connected component. 

In the following, we introduce some basic concepts and results of CAD. The reader is 



referred to Collins (1975), Hong (19901, McCallum (1988 1998), Brown (2001) and Xiao 



(2009) for a detailed discussion on the properties of CAD and Open CAD. 



Definition 16. (Collins 1975 McCallum 19881 An n-variate polynomial f{X n -\, x n ) 
over the reals is said to be delineable on a subset S (usually connected) of M™ -1 if (1) 
the portion of the real variety of / that lies in the cylinder SxR over S consists of the 
union of the graphs of some k > continuous functions B\ < ■ ■ ■ < 8^ from S to K; and 
(2) there exist integers m\ , . . . , mk > 1 such that for every a E S, the multiplicity of the 
root 9i{a) of f{a,x n ) (considered as a polynomial in x n alone) is rrij. 
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Definition 17. (Collins, 1975 McCallum 1988) In the above definition, the Oi are called 



the real root functions of / on S, the graphs of the &i are called the /-sections over S, 
and the regions between successive /-sections are called /-sectors. 



Theorem 18. (McCallum, 1988\ 1998) Let f(X n ,x n +i) be a polynomial in R[X n , x n +i] 
of positive degree and discrim(/, x n +i) is a nonzero polynomial. Let S be a connected 
submanifold of R" on which f is degree-invariant and does not vanish identically, and 
in which discrim(/, x n+ i) is order-invariant. Then f is analytic delineable on S and is 
order-invariant in each f -section over S. 

Based on this theorem, McCallum proposed the projection operator MCproj, which 
consists of the discriminant of / and all coefficients of /. 



Theorem 19. (Brown 2001) Let f(X n , x n+ i) be a (n + l)-variate polynomial of positive 
degree m in the variable x n +i with discrim(/, x n+ i) ^ 0. Let S be a connected submanifold 
of M. n in which discrim(/, x n +i) is order-invariant, the leading coefficient of f is sign- 
invariant, and such that f vanishes identically at no point in S. f is degree-invariant on 
S. 

Based on this theorem, Brown obtained a reduced McCallum projection in which only 
leading coefficients, discriminants and resultants appear. The Brown projection operator 
is defined as follows. 



Definition 20. (Brown 2001) Given a polynomial / € K[A"„ + i] with level n + 1, the 



Brown projection operator for / is 

u ( sqrfrcc ( / ) ) 

Bproj(/,x n+ i) = Res(sqrfree(/), ,x n+ i). 

If L is a polynomial set and the level of any polynomial in L is n + 1, then 

d ( scjrf ree ( /* ) ) 

Bproj(L,x I1+ i) = (J {Res(sqrfree(/), — ,x n+1 )} 



feL 



dx n+1 



|J |J {Res(sqrfree(/),sqrfree(5),a;„ + i)}. 

f,g£L,f^g 



Algorithm 1. Bprojection ( |Brown| |2001[ ) 
Input: A polynomial f{X n ) S M[X„]. 
Output: A projection factor set F. 



F := {/(*„)}; 

for i from n downto 2 do 

F := FU{Bproj(F l ,Xi)}; (F l is the set of polynomials in F with level i). 
end for 
return F 

Open CAD is a modified CAD construction algorithm, which was named in Rong 



Xiao's Ph.D. thesis (Xiao 2009). In fact, Open CAD is similar to the Generic Cylindri- 



cal Algebraic Decomposition (GCAD) proposed by Strzebohski (2000) and was used in 



G 



DISCOVERER (Xia 2000) for real root classification. For convenience, we describe the 
framework of the Open CAD here. 

For a polynomial f(X n ) £ M[X n ], an Open CAD defined by f(X n ) is a set of sample 
points in R™ obtained through the following three phases: (1) Projection. Use the Brown 
projection operator (Algorithm 1) on f(X n ); (2) Base. Choose one point in each of the 
open intervals defined by the real roots of F 1 (see Algorithm 1); (3) Lifting. Substitute 
each sample point of R I_1 for in F l and then, by the same method as Base phase, 

choose sample points for F l {xi). 



The Successive Resultant Method 



The Successive Resultant Method (SRes) is introduced in (Yang 20011 without a 



proof. The method can be used for solving Problem 2 of this paper, i.e., problem of 
global optimization. 

For a polynomial f(X n ) in Problem 2, the SRes method first applies Brown's projec- 
tion recursively on polynomial f{X n ) — K to get a polynomial g{K). Suppose g(K) has m 
distinct real roots ki(l < i < to). Then computes to + 1 rational numbers Pi(0 < i < m) 
such that ki £ (pi_i,pi). Finally, substitutes each pi for K in f(X n ) — K to check if 
f(X n ) — pi > holds for all X n € R. If pj is the first such that f(X n ) — pj > does not 
hold for all X n £ R, then kj is the infimum (fco = -co). To check if f(X n ) — pi > holds 
for all X n £ R, the SRes method applies Brown's projection on f(X n ) — pi and choose 
sample points only from the highest dimensional cells in the lifting phase. 

The SRes method is formally described as Algorithm 2 and we will prove its correctness 
in the rest parts of the section. 



Algorithm 2. SRes (Successive Resultant Method, [Yang] \2QQl\ ; |Yang et al.| ( |2008| ) 
Input: A squarefree polynomial / £ R[X„]. 

Output: The supremum of k £ R, such that VX n £ R™, f(X n ) > fc. If there doesn't 

exist such fc, then returns — oo. 
l: F :— Bprojection(f — k) (F l is the set of polynomials in F with level i. Here F % has 

no more than one polynomial, we denote this polynomial by 
2: Co :=an Open CAD of R defined by F (Suppose Co = \JiL Pi, Pi & {h, h+i), where 

ki (1 < i < to) are the real roots of F and k = — oo, k m+ i = +oo.); 
3: for I from to to do 
4: for i from 1 to n do 

5: Cu~ an Open CAD of R l defined by [j] 

6: end for 

7: if 3X n £ C ln , s.t. f(X n ,pi) < then 
8: return ki 

9: end if 
10: end for 



=1 Fj and k 



Pi, 



Remark 21. If g(X n ) > for all X n £ R™, Algorithm 2 can also be applied to compute 
inf{ g(x") \Xn € R™}. We just need to replace F := Bprojection(f — k) of step 1 by 
F := Bprojection(f — kg). The proof of the correctness is the same. 
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We first prove the following lemma, which can be inferred from the results of ( Brown 



2001) and (McCallum 1998), namely / is delineable over the maximal connected regions 
defined by Bproj(/, x n ) 7^ 0. Since the proof is inspiring, we will present it here. 



Lemma 22. 



( McCallum\ \ 1998\ - \Brown\ \2001) Let F iy 

l 



_i be as in Algorithm 2. Let U 
be a connected component of -Fi-i 7^ in K I_1 and y}{^) < yf(j) < ■ ■ ■ < y™{l) be all 
real roots 0/^(7, Xi) = for any given 7 £ U. Then for all a, (3 £ U , a x (yf (a), yl (a)) 
and (3 x Q/f ~ y{{P)){j = 2, 3, . . . , m) are in the same connected component of Fi 7^ 
in W . Therefore, for any two points in U, their corresponding sample points obtained 
through the Open CAD lifting phase are in the same connected component of F n 7^ in 



Proof. For a £ U, let e = min | yi(a) — yi-i(a) |, by Lemma 15, 3d > 0, such that 

2<i<m 



Va' £ B{a,8), max | y t (a) - yda') |< §. 

l<i<m 

Consider segment (a, ^-i(<*)+fe(°0 ) ^-1 (»')+» ("') ) where a / £ £(a,S). For 

any point (a",j/) on the segment, we have 

I ¥-».(«") I 

= I (y s («") - *(«)) + („.(«) - ^- l(Q) 2 + ^ (a) ) + ( Ifc^tltfW _ y) 1 

> I ».(«) - ^toW 1 _ 1 y s{a ")-y s{ a) | - | j^t^M . y | 

~ 2 6 1 2 2 1 

see 

~2 ~ 6 ~ 6 
>0 

So the points satisfying Fi = are not on the segment. 

Therefore, for any points r 1 £ (y j _ 1 (a),y j (a)),r2 £ (y j ^ 1 (a'),y j (a')), (a' £ B a (5)), 
the points satisfying Fi — are not on the broken line (a,ri) — > (a, o" 1 " 2 ' 3 ) — ► 

( a ' /'- l( °'^ (n V (aV2). 

Hence we know that for any a € £7, there exists $ > such that for any point 
a' £ B a (S) and 2 < s < m, a x y s (a)) and a' x (j/ a _i(a'), y s (a')) are in 

the same connected component of Fi 7^ in W. For all a,/3 <G U, there exists a path 
7 : [0, 1] — > U that connects a and /3. Due to the compactness of the path, there are finitely 
many open sets B at (5 t ) covering 7QO, 1]) with a t £ 7([0, 1]), Va' £ B at (5 t ),2 < j < 
m,ax (jjj_i(a),yj(a)) and a' x (j/j_i(a'), yj(a')) are in the same connected component 
of Fi 7^ 0. Since the union of these open sets are connected, the lemma is proved. □ 



Remark 23. By the above Lemma, in Algorithm 2, for any two points pi, p[ £ (ki, &z+i), 
their corresponding sample points obtained through the Open CAD lifting phase are in 
the same connected component of F n 7^ in R n+1 . Since at least one sample point can 
be taken from every highest dimensional cell via the Open CAD lifting phase, the set 
of the corresponding sample points of pi obtained through the Open CAD lifting phase, 
Ci n in Algorithm 2, contains at least one point from every connected component U of 
F(X n ,k) ^ 0, in which UC]M. n x {h-\,h) ^ 
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Theorem 24. The Successive Resultant Method is correct. 

Proof. Let notations be as in Algorithm 2. If there exists a A/ € such that 

F(X n , k') > for all X n £ MP, then by Lemma 22, for any k £ (k u k i+1 ), F(X n , k) > 
for all X„ G R" (since their corresponding sample points obtained through the Open 
CAD lifting phase are in the same connected component of F n ^ in W 1 ). Therefore, 
for any k € [fc i; fc i+1 ], F(X n ,k) > for all X n £ R™. The global optimum k will be 
found by checking whether VA„ £ IR" - , F(X n ,pi) > holds where ^ is the sample point 
of (k i7 k i+ i). Since Algorithm 2 ensures that at least one point is chosen from every 
connected component of F n (X n ,pi) ^ in R n , the theorem is proved. □ 

4. Solving Problem 1 via simplified CAD projection 

To solve Problem 1, i.e., to prove or disprove f(X n ) > 0, one may take the following 
two steps: choose at least one point from each connected open cell defined by f(X n ) in 
M™ and then check whether or not f(X n ) > at all those sample points. A CAD based 
method for choosing those sample points need to do projection by, for example, the 
Brown operator. So, we need to compute discriminant, discriminant of discriminant, and 
so on. Now, suppose f(X n ) is semi-definite, by Lemma 13, discrim(/, x n ) and lc(/, x n ) 
are both semi-definite. By using Lemma 13 again and again, we know that all polynomials 
appearing in the Brown projection are semi-definite. 

Based on this observation, in this section, we will improve Brown's projection for solv- 
ing Problem 1 from the perspective of maintaining scmi-dcfinitcncss. A new projection 
operator Nproj is proposed, which avoids some resultants calculation and thus reduces 
the scale of projection. 

4-1. Notations 

Definition 25. Suppose h £ M[X n ] can be factorized in M[A„] as: 

h = i^-Hl*- 1 . . . i^h^h^ . . . h m 21 -, 

where hi(i — 1, 2, . . . , m) and = 1,2, ... ,t) are pairwise different irreducible polyno- 
mials in K[X„]. We define sqrfree 1; sqrfree 2 as follows. 

sqrfree 1 (ft) = {k,i = 1,2,...,*}, 
sqrfrcc 2 (/i) = {hi,i = 1,2,..., m}. 

Suppose / £ R[X n ] is (positive/negative) semi-definite. To make the following de- 
scription simpler, we assume / is monic w.r.t. the variable x n . Let h = discrim(/, x n ). As 
discussed at the beginning of this section, we have to choose sample points from every 
open connected component defined by sqrfrcc(/i) = sqrfrcc 1 (/i) (J sqrfree 2 (ft.). By Lemma 
13, discrim(/, x n ) is semi-definite. We will prove that the polynomials in sqrfree 1 (/i) are 
also semi-definite. It will also be shown that for any semi-definite polynomial g £ M.[X n ], 
]R™\Zero(gf) is connected. This property guarantees that the open connected components 
defined by sqrfree(/i) = sqrfree 1 (/i) \J sqrfree 2 (/i) are almost the same as the open con- 
nected components defined by sqrfree 2 (/i) (may differ with the real roots of sqrfree 1 (/i)). 
We call a point is "good" if it is not a zero of any polynomial in sqrfree 1 (ft,). There- 
fore, in order to prove / > for all X n £ R™, we only need to choose "good" sample 
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points from every open connected component defined by sqrfree 2 (/i) and prove that the 
polynomials in sqrfree 1 (ft.) are semi-definite. The latter can be solved by recursion. Thus 
in the CAD projection phase, the resultants of those polynomial pairs of which one is 
from sqrfree 1 (/i) and the other from sqrfree 2 (/i), as well as the resultants among the 
polynomials in sqrfree 1 (/i) are unnecessary. 

Having this observation in mind, we propose a new projection operator Nproj, the 
scale of which is no larger than that of Brown's. Based on the new operator, we obtain 
a new algorithm Proineq to prove or disprove inequalities. 

We first introduce the principal and secondary parts of the new projection operator 
Nproj . 

Definition 26. Given a polynomial / e R[A„ +1 ] with level n + 1. We define 

cs(/,ar n+ i) = sqrfree 2 (lc(/,ar n+ i)),ms(/,a; n+ i) = sqrfree 2 (discrim(/, x n+l )). 
The secondary and principal parts of the new projection are defined as 

N P ro j a (/> x n+i) = sqrfree! (lc(/, x n+1 ) ) (J sqrfreej (discrim(/, x n+1 )) , 

Nproj p (/,a; n+ i) = { J[ g}. 

a 6 cs(/, x n+1 ) ms(/, x n+ i), 
g £ nproj s (/,a;„ + i) 

While L is a set of polynomials with level n + 1, we define 

Nproj s (L,a;„ +:L ) = U /eL sqrfree 1 (lc(/,x„ +1 ))(Ju /eL sqrfree 1 (discrim(/,x„ +1 )), 

Nproj p (L,a;„ + i) = (J { J| g}. 

f eL S effl(/,l„ +1 )Um.(/,i„ +1 ), 
g & Nproj 3 (L, x n+1 ) 

Now the algorithm Nproj is clear and is formally described here. 

Algorithm 3. Nproj 

Input: A polynomial f(X n ) € 

Output: A projection factor set (Li,i 2 ). 

l: L 2 := sqrfree 1 (/); 

2: L x := {}; 

3: for i from n downto 2 do 

4: L\ := L\ |J Nproj p (L 2 , Xi) [j Bproj(L^, Xj); (Notice that L % is the set of polyno- 
mials in L with level i.) 
5: L 2 :=L 2 \JHpToi 3 {L} l ,x i y, 
6: end for 
7: return (Li, L 2 ) 

4-2. The correctness of Algorithm Proineq 

Theorem 27. Let f and g be coprime in M.[X n ]. For any connected open set U in M. n , 
the open set V = ?7\Zero(/, g) is also connected. 
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Algorithm 4. Proineq 

Input: An irreducible polynomial / £ R[A„]. 
Output: Whether or not VA £ R", f(X n ) > 0. 

1: (Li,L 2 ) :=Nproj(/) 

2: for i from 1 to n do 

3: Ci~ An Open CAD of W defined by U} =1 Z/[ (J U} =1 ^; 
4: if 3A/, A? £ Ci and 5 £ L\, s.t., g(Xl)g(Xf) < then 
5: return false 

6: end if 

7: if 3A* £ C n , s.t., /(A*) < then 

8: return false 

9: end if 
10: end for 
11: return true 

Proof. For any two points a, [3 in V, we only need to prove that there exists a path 
r(t) : [0, 1] — > V such that r(0) = a, r(l) = /3. Choose a path ru that connects a and f3 
in 17. Notice that 17 is an open set, so for any A„ £ ru, there exists 5x n > such that 
17 D -Bx„(<5x„)- Since is compact and U-^^n(^n) i s an open covering of ru, there 

m 

exists an in £ N, such that (J B X k(5xk) D ru and a £ Bxj(ixi), /? £ Bx™(5x™)- Now 

A— 1 

we only need to prove that B X k(5 X h)\Zevo(f , g) is connected. 

Choose a, & £ 5^ (<^x*)\Zero(/, 5). There exists an affine coordinate transformation T 

such that T{B X k{8 X k)) = B 0n (l) and T(a)T(b) and (0„_i,l) are parallel. Thus the first 
n— 1 coordinates of T(a) and T(6) are the same. Let T(a) = a'), T(b) = {Y n _x, b'). 

Without loss of generality, we assume that a' > b' . 

In the new coordinate, / and g become T{f) and T(g), respectively. i?o„(l) is an 
open set and T(a),T(b) £ Zevo(T(f),T(g)), so there exists r > such that the cylin- 
der By^ir) x [b',a'} C S 0b (1), B T(a) {r) f| Zero(T(/), T(g)) = and B T(b) (r)f] Zero 
(T(f),T(g)) = 0. By Lemma 9, T(/) and T(g) are coprimc in R[A„]. So by Lemma 10, 
there exists X' n _ 1 £ By n _i(r), such that for any x n £ R, (A"4_i,a;„) ^ Zero (T(f),T(g)). 
Thus the broken line T(a) — >• (A^-^a') — > (A^-^fc') — > T(6) is a path that connects 
T(a) and T(6) in B 0n (l)\Zero(T(/),T(p)). The theorem is proved. □ 

Proposition 28. Let U C M™ &e a connected open set, f, g £ R[A„], (/, ff) = 1 and 

/or aZZ A„ £ 17, f(X n )g(X n ) > 0. CTfcen f(X n ) > 0,g(X n ) > for all X n £ 17 or 
f(X n ) < 0,g(X n ) < for all X n £ U. 

Proof. If not, there exist A*, A 2 £ 17, such that /(A*) < 0, g(A'^) < and /(A 2 ) > 0, 
g{Xn) > 0. By Theorem 27, U\Zero(f,g) is connected. So we can choose a path 7 that 
connects A* with A 2 and 7 1") Zcro(/, g) = 0. Consider the sign of / + g on 7. Since the 
sign is different at A* and A 2 , by Mean Value Theorem we know there exists A^ on 7 
such that f(X^) +g(X%) — 0. From the condition we know that f(X n )g(X n ) > 0, hence 
X^ £ Zero(/, g), which contradicts the choice of 7. □ 
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Proposition 29. Given a non-zero squarefree polynomial f G R[X n ] and a connected 
open set U ofR n , if f(X n ) > for all X n G U, then £/\Zero(/) is also a connected open 
set. 

Proof. Since / is non-zero, we only need to consider the case that the level of / is non- 
zero. Let i > be the level of / and consider / as a polynomial of Xi. Because f{X n ) > 
on U, we know Zero(/) f| U = Zero(/, f' x .) f| U. 

If / is irreducible in K[X n ], / and f are coprime in RLY„]. Thus C/\Zero(/, f' x .) is 
connected by Theorem 27. So C/\Zero(/) is a connected open set. 

If / is reducible in K[X„], let / = Yit=i /*> m which & h ft(t = 1, • • • ,j) are irreducible 
polynomials in then C/\Zero(/) = U\ U*=i Zero(/ t ) is a connected open set. The 

proposition is proved. □ 

Theorem 30. Given a positive integer n > 2. Let f G KLY„] be a non-zero squarefree 
polynomial and U be a connected component o/Nproj p (/, x n ) ^ in W 1 ^ 1 . If the polyno- 
mials in Npro j s (f,x n ) are semi- definite on U , then for a, j3 G U\ UheNproj (fx ) Zero(/i), 
t/ie corresponding points in adjacent real roots of f(a,x n ) = anrf f(f3,x n ) = are m 
t/ie same connected component of f ^ 0. 

Proof. By Proposition 29, U\ U/ieNproj (/ k ) Zero(/i) is also a connected open set. So a, 
j3 are two points in the same connected component of Bproj(/, x n ) ^ 0. By Lemma 22, 
the corresponding points in adjacent real roots of f(a,x n ) = and f(/3,x n ) = are in 
the same connected component of / ^ 0. □ 

Theorem 31. Given a positive integer n>2. Let f G KLY„] be a squarefree polynomial 
with level n and U be a connected open set o/Nproj p (/, x n ) ^ in The necessary 

and sufficient condition for f(X n ) to be semi-definite on U x K is the following two 
conditions hold. 

(1) The polynomials in Nproj s (/, x n ) are semi-definite on U. 

(2) For an arbitrary point a G C^\Uft,eMproj (/x„) Zero(ft), f(a,x n ) is semi-definite on M. 

Proof. =>: By Lemma 13, discrim(/, x n ) is semi-definite on U. Thus by Proposition 
28, the polynomials in Nproj s (/, x n ) are semi-definite on U. It is obvious that f(a,x n ) 
is semi-definite on R. 

<;^: If the polynomials in Nproj s (/, x n ) are semi-definite on U, by Theorem 30, 
^\ U/ieNproj (/ x ) Zero(ft) is a connected open set. From that f(a, x n ) is semi-definite on 
R, we know that f(X n ) is semi-definite on U x R. □ 

The following theorem is an easy corollary of the above theorem. 

Theorem 32. Given a squarefree polynomial f G with level n, the necessary and 

sufficient condition for f{X n ) to be semi-definite on R n is the following two conditions 
hold. 

(1) The polynomials in Nproj si (/, x n ) are semi-definite on R™ -1 . 

(2) If cti(i = 1,2) are two arbitrary points from any two different connected components 
of Nproj p (/, x n ) in and oti is not a zero of any polynomial in Nproj g (/, x n ), 
then f(a>i, x n ) > (i = 1,2) for all x n G M or /(a,, x n ) < (i = 1, 2) /or aiZ i„ G R. 
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Similar to the CAD lifting phase described in ( Arnon et al. 1984) and (Brown 2001 1 



Algorithm 4 proves (or disproves) the inequality f(X n ) > by testing the sample points 
computed by Open CAD. From Theorem 32, we know that Algorithm 4 is correct. 



5. Solving Problems 2 and 3 

This section mainly considers Problem 2 and Problem 3 put forward in the intro- 
duction. Without loss of generality, we only need to consider the case that the given 
polynomial is squarefree. 

We first propose an algorithm, named Findk, for solving Problem 3, which is based 
on our new projection operator Npro j . For a given f(X n , k), by Npro j (/), we can obtain 
a polynomial P(k) in k. Assume fcj(l < i < m) are the real roots of P(k). We only need 
to check whether or not f(X n , k) > holds at those fc^s and on the m + 1 open intervals. 
The algorithm is formally described as Algorithm 5 and its correctness will be proved in 
Theorem 34. 

Lemma 33. In Algorithm 5, 

(1) For any polynomial fi(Xi,k) in L\, Cu contains at least one point from every con- 
nected component U of ^ 0, in which U f]W x (fc;_i,fej) ^ 0; 

(2) When < i < n, for any polynomial gi(Xi,k) in L\, if pi £ GKf, then gi(Xi,k) is 
semi-definite on W x (/c ; _i, k t ). Herein, GKf — FKf U FK_f. 

Proof. (1) It can be easily deduced from Remark 23. 

(2) We prove it by induction on i. When i = 0, the conclusion is obvious. When i = 1, 
by Theorem 31, it is also true. Assume the conclusion is true when i — j — l(j > 2). For 
any polynomial gj(Xj,k) in L 2 , notice that Nproj s (gj) C L 2 , Nproj p (<?j) C L\ .If 
Pi 6 GK gj , from the process of the algorithm, we know that pi £ GK g for any polynomial 
g in Nproj s {gj). So by the assumption of induction, we know that every polynomial in 
Nproj^cfy) is semi-definite on K^ -1 x (fc;_i,fc;). By (1) and Theorem 31, gj(Xj,k) is 
semi-definite on K J x (fc;_i,fc;). That finishes the induction. □ 

Theorem 34. The output of Algorithm 5, FKf, is {k\VX n E R n ,f(X n ,k) > 0}. 

Proof. Denote {k\\JX n e W\ f(X n ,k) > 0} by K f . 

We first prove that FKf C Kf. It is sufficient to prove that if pi € GKf, then / is 
semi-definite on R™ x (fc;_i, k[). Notice that / € L 2 , so this follows from Lemma 33(2). 

We then prove that Kf C FKf. It is sufficient to prove that if there exists k' € 
(h-uh) such that VX„ € R", f(X n , k') > 0, then (fcj_i,fcj) € FKf. 

It is obviously true when n = 1. When n > 2, for any polynomial g n —i(X n —i,k) 
in sqrfree!^-!^-!,^)) C Nproj S (/(X„, k'),x n ). If f(X n ,k') is squarefree, 

because VX n e R n ,f(X n ,k') > 0, by Theorem 32, <7 n _i(Jf„_i, fc') is semi-definite on 
1R 71-1 . When f(X n , k') is not squarefree, g„_i(X n _i, fc') = and it is also semi-definite. 
So g n -i(X n —i, k 1 ) is semi-definite on In a similar way, we know that for any 

1 < j < n — 1 and any polynomial gj(Xj, k) in L J 2 , gj(Xj, k) is semi-definite on M J x fc'. 

If the algorithm is incorrect, then there exists a minimal i and gi(Xi,k) £ L 2 , such 
that pi GK gi . If such gi is not unique, then choose any one. If i = 1, it contradicts 
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Algorithm 5. Findk 

Input: A squarefree polynomial f(X n , k) e k}. 
Output: The set FK f = {k\iX n e K n , f(X n , k) > 0}. 
l: FKf := {}; 

(Li,L 2 ) :=Nproj(/); 

Co := a CAD of R defined by L?U L °; (Suppose C = U^fe} {fcj where 
< i < m) are the real roots of L\ U = — k m+ \ = +oo, pi e fcj+i). 

is the set of polynomials in L with level i.) 
for I from to m do 

Choose sample point pi of (fc/, 
w := 1; 

for i from 1 to n do 

C (l :=an Open CAD of W defined by U* =1 Lj (J U) =1 L 3 2 ,k = Pl ; 
if e C M , 5 e L\, s.t. g{X} ^giXlm) < then 

i> := 0; 
break 
end if 

if i = n and 3X n e Cj n , s.t. f(X n ,pi) < then 

i> := 0; 
end if 
end for 
if f = 1 then 

FK f :=FK f \J(k l ,k l+1 ); 
end if 
end for 

for i from 1 to m do 

if Proineq(f (X n , fcj)) =true then 

f^ /: =f^u{M; 

end if 
end for 
return FKf 



Theorem 31. Hence i > 2. From the definition of i we know that G Gifg for any 
polynomial g in Nproj s (gj). By Lemma 33(2), every polynomial of Nproj s ((^ i ) is semi- 
definite on x (ki-i,ki). By Lemma 33(1), for any polynomial in L\~ , C;, 
contains at least one point from every connected component U of 7^ in which 
[/f|K i_1 x (fcj_i,fcj) ^ 0. Thus, by Theorem 30, the indefiniteness of gi{X u k) on the 
sample points of W x (fc;_i,fc;) will lead to the indefiniteness of gi(Xi,k'), which is a 
contradiction. The proof is completed. □ 

By modifying the algorithm Findk, we get the algorithm Findmax which can solve the 
global optimum problem (Problem 2). 

6. Examples 

In this section, we will illustrate the process and efficiency of algorithms Findmax, 
Findk and Proineq with several non-trivial examples. 
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Algorithm 6. Findmax 
Input: A squarefree polynomial / e 
Output: iel such that k = mix gr™ f{X n ). 
L:={} 

(Li,L 2 ) :=Nproj(/(X n )-fc) 

Co :=a CAD of K defined by L° |J L°. (C = U^fe} (J U^Jfc;}, k = -oo,k m+1 = 

+00, pi e (fci,fc i+ i)) 

for ^ from to m do 

Choose sample point p; of (fc;, 
v:=l; 

for i from 1 to n do 

C (l :=an Open CAD of W defined by U} =1 ij (J Uj =1 L 3 2 and fc = p ; ; 
if 3Ai,A t 2 e C u ,g € L\, s.t. g{X}, k)g(Xf, k) < then 
v := 0; 
break 
end if 

if i = n and 3X n <G C7 n , s.t. f(X n ) — p t < then 

i> := 0; 
end if 
end for 
if v = then 
return fc ; 
end if 
end for 
return k m+1 

Remark 35. For f,g e K[A„], VA„ e E™, e/(X„) > 0, algorithm Findmax can also 
be applied to compute inf{ ^^"| \X n G E"}. We just need to replace (Li,L 2 ) := 
Nproj(/(A„) - fc) of step 2 by (Li,L 2 ) := Nproj(/(X„) - fc<?(X„)). 



Algorithm Findmax will be compared with the algorithm SRes. The program Proineq 
we implemented using Maple will be compared with the function PartialCylindricalAl- 
gebraicDecomposition of RcgularChains package in Maplel5. 

All computations were performed on a laptop with Inter Core2 2.10GHz CPU and 
2.00GB RAM. 

Since our main contribution is the new projection operator Nproj, for the following 
example, we will mainly focus on how the projection phase of algorithm Findmax differs 
from that of algorithm SRes. 

Example 6.1. Compute inf^^^gR G(x, y, z), where 

c ^ (x 2 -x + l)(y 2 -y + l)(z 2 -z + l) 
(xyz) 2 — xyz + 1 

Since for any x, y, z e E, 

(xyz) 2 — xyz + 1 > 0, 
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so this problem can be solved either by algorithm Findmaxor algorithm SRes. We will 
apply algorithm Findmax and algorithm SRes to solve this problem, respectively. 

1. Algorithm Findmax. 

Projection phase. 

Notice that in the projection phase of Findmax, we call algorithm Npro j . 
For the first step, L 2 := {/}, L\ := {}, where 

/ := (x 2 -x+ l)(y 2 - y + l)(z 2 - z + 1) - k((xyz) 2 - xyz + 1). 

Now we step into the loop of algorithm Npro j . 

When i := 3, we regard / as a polynomial in x. Then lc(/, x) = /i, discrim(/, x) = f 2 . 
So L\ := {}, L 2 := {/i,/ 2 }, where 

h :=(y 2 -y + l)(z 2 -2 + l)-%V, 

h ■= - 3(1 -2z-^k-2y + 3y 2 + 3z 2 - ^kyz + 9y 2 z 2 

- 6y 2 z - 6yz 2 + Ayz - 2ky 2 z 2 + z 4 y 4 - 2y 4 z 3 
+ 3y 4 z 2 - 2y 3 z 4 + Ay 3 z 3 - 6y 3 z 2 + 3y 2 z 4 

- 6y 2 z 3 - 2y 4 z + Ay 3 z - 2yz 4 + Ayz 3 + ^kz 
+ \ky - ^ky 2 - ^kz 2 + k 2 y 2 z 2 

+ g«y z + T^kyz - -kz y 

+ 4 ky 4 z 3 - % 4 z 2 + hy 3 z 4 
3 3 y 3 y 

- -ky 3 z 3 + -ky 3 z 2 - -ky 2 z 4 
3 3 y 3 y 

+ 2 -ky 2 z 3 + ^ky 3 z + ^kyz 3 + y 4 
+ z 4 - 2z 3 - 2y 3 ). 

When i := 2, we regard /i, f 2 as polynomials in y. Then lc(/i, y) = /n, discrim(/i , y) = 
A2/13, lc(/2,J/) = /21/22, discrim(/ 2 ,y) = 2304/f 1 / 22 /2 3 /24- Notice that / 13 = / 2 i, so 
:= {/24L ^2 : = {/u,/i2,/i3,/22,/23}, where 

/u :=-(fcz 2 -z 2 + z-l), 
/12 := - 3z 2 + 4fcz 2 + 3z - 3, 
/is :=^ 2 -2 + 1, 
/21 :=^ 2 -2 + 1, 

/ 22 :=9z 4 fc 3 - 24z 6 fc 2 - 21z 5 fc 2 - 3z 4 k 2 - 2\z 3 k 2 

- 2Az 2 k 2 + 16fcz 8 - Akz 7 + 16fcz 6 + 32fcz 5 - 8kz 4 

+ 32/cz 3 + 16fcz 2 - Akz + 16fc - 9z 8 + 9z 7 - 9z 6 - 18z 5 

+ 18z 4 - 18z 3 - 9z 2 + 9z - 9, 

/ 23 := - 1 + 2z - 3z 2 + 2z 3 - z 4 + fcz 2 , 

/24 :=fc- 
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When i := 1, we regard /n, /12, /13, /22 and /23 as polynomials in z. Then lc(/n,z) = 
1 — fc, discrim(/n, z) = — 3 + 4fc; lc(/i2,^) = — 3 + 4fc, discrim(/i 2 , 2:) = 3(— 9 + 16k); 
lc(/ 13 ,z) = 1, discrim(/ 13 ,z) = -3; lc(/ 22 ,*) = -9 + 16fc, discrim(/ 22 , z) = 531441(-1 + 
fc)(fc - 9)(3fc - 2) 2 (fc 2 + 6fc - 3) 2 (256fc - 279) 6 fc 14 ; lc(/ 23 ,z) = -1, discrim(/ 23 , z) = 
16fc 2 (-l + fc)(fc-9). So 

L? = {3fc - 2, fc 2 + 6fc - 3, 256fc - 279, fc}, L% = {1 - fc, -3 + 4fc, -9 + 16fc, fc - 9}. 

Thus, the projection phase is finished. We get a polynomial in k defined by L\ (J L 2 , 

g := (3fc - 2)(fc 2 + 6fc - 3)(256fc - 279)fc(l - fc)(-3 + 4fc)(-9 + 16fc)(fc - 9). 

By real root isolation, we get 10 sample points defined by (J L 2 , 

r r 1 7 n r r 1 , ,-17-, r 1 9 -, r 23 -, r 7 -, r 33 -, r 73 -1 r 33 11 

[[^],[-2],[3],y,y,y,y,y,y,[ T ]]. 

After testing these points, we know that the maximum k is the real root of k 2 + 6k — 3, 
which is between | and |. 

2. Successive resultant method 

Projection phase. 

Notice that in the projection phase of the Successive Resultant Method, we call algo- 
rithm 1 (Bprojection). 

In the first step, F := {/}, where 

/ = (x 2 -x + l)(y 2 - y+ 1)0 2 - z + 1) - k((xyz) 2 - xyz + 1). 

Now we step into the loop of algorithm Bprojection. 
When i := 3, we regard / as a polynomial in x. Then 

Res(f,f,x) = f 1 f 2 , 

F 2 := sqrfree(Res(/, /', a;)) = fif 2 . 

When i := 2, we regard / as a polynomial in y. Then 

Res(F 2 ,F^y) = 2304fc 12 z 8 /n/ 1 2 2 / 1 8 3 /2 2 / 2 7 3 , 

Fx := sqrfree(Res(F 2 ,i^,?/)) = kzf n fi 2 f 13 f 22 f 23 . 

When i := 1, we regard / as a polynomial in z. Then 

Res(Fi,i^,z) = - 1156831381426176fc 125 (-9 + 16fc) 4 (-3 + 4fc) 2 
(fc 4 - 294fc 3 + 1425fc 2 - 2277fc + 1089) 2 (-1 + fc) 5 
(614656fc 4 - 4409856fc 3 + 11013408fc 2 - 11477376fc + 4021893) 2 
(3ft - 2) 2 (fc 2 +6k- 3) 2 (256fc - 279) 6 (fc - 9) 14 (4fc - 9) 4 , 

F :=sqrfree(Res( J Fi, F{, z)) = g(k 4 - 294fc 3 + 1425fc 2 - 2277fc + 1089)(4fc - 9) 
(614656fc 4 - 4409856fc 3 + 11013408fc 2 - 11477376fc + 4021893). 

The maximum fc is a real root of Fq. By real root isolation, we get 15 sample points of 
fc defined by F , 

rr -17. , „. .7. r 33 n r 39 n r 87 n r 47 n r 25 n r 15. r 33 n r 25 n ,11, r 23 n , „„. 4025.. 
[[^-], [-2], [gj], y , [&]> ^ [ 32 ] ' [ 16 ] ' [ T ] ' [ T ] ' [136] ' [ — «■ 
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After testing these points, we know that the maximum k is the real root of k 2 + 6k — 3, 
which is between | and |. 

Obviously, the scale of projection with the new projection is smaller. The polynomial 
in k calculated through the successive resultant method has three extraneous factors. 



Example 6.2. Let X± € K , prove or disprove 

4 



i=l 



3017 
1000 



^xfzj+i > 0, 



where X5 = Xi. 

Notice that /(l, 1, 1, 1) = ||§ > 0, it is a non-trivial example. When we use Proineq 
to prove or disprove / > 0, after 0.921 seconds, Proineq returns false. If we use Partial- 
CylindricalAlgebraicDecomposition, it needs 2.187 seconds. 

Actually, the maximum value of k satisfying the following inequality, 

4 4 

(VX 4 e K 4 ) f(X 4 , k) = (J2 x 2 ) 2 - k x\xi+x > 0, 



8 j 40 



= 1 
3016 



is —7^ + v 5, which is between |^ and 1( 



Example 6.3. (Han 2011) Let I„ € R", then 



F(X n ,n)=l[( 



x 2 + n-l)-n n - 2 



£>) 2 >0. 



When n = 3,4,5,6, we compared PartialCylindricalAlgebraicDecomposition(PCA£>) 
with Proineq. 
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PCAD 




Proineq 






Time 


Number of 


Time 


Number of 




(seconds) 


sample points 


(seconds) 


sample points 


3 


0.078 


4 


0.063 


3 


4 


0.250 


10 


0.422 


3 


5 


2.282 


36 


0.875 


5 


6 


>3000 


Unknown 


4.188 


15 




Example 6.4. Let X n e 


R n , n = 3,4,5, 


then 





where x n +i — x\. 

We notice that the inequality is homogenous, thus we can assume that x n = 1. Now 
we compare PartialCylindricalAlgebraicDecomposition(PCAZ)) with Proineq. 
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n 


PCAD 




Proineq 






Time 


Number of 


Time 


Number of 




(seconds) 


sample points 


(seconds) 


sample points 


6 


U.Uoz 





U.U4 I 


r. 



4 


0.469 


2 


0.171 


2 


5 


>3000 


Unknown 


244.188 


20 




Example 6.5. Let 


yi>0,i = l,2,...,8 


, prove that 








2.2.2 
" + Vi+2 < Vi+3 ~ 


2.2,2 
Ui+4 + Vi+5 + Vi+6 


- y 2 i+ 7)) > °> 



i=l 



where y j+s = yj for j = 1,2,..., 7. 

It takes 27.218 seconds for Proineq to prove that F > after checking 6720 sample 
points. PCAD cannot obtain result on this example within one hour. 

The above three examples demonstrate that in terms of proving inequalities, algorithm 
Proineq is efficient. Further improvements on the projection and lifting phase are our 
ongoing work. 
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